Phases and collective modes of hardcore Bose-Fermi mixture in an optical lattice 
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We obtain the phase diagram of a Bose-Fermi mixture of hardcore spinless Bosons and spin- 
polarized Fermions with nearest neighbor intra-species interaction and on-site inter-species repulsion 
in an optical lattice at half-filling using a slave-boson mean-field theory. We show that such a system 
can have four possible phases which are a) supersolid Bosons coexisting with Fermions in the Mott 
state, b) Mott state of Bosons coexisting with Fermions in a metallic or charge-density wave state, 
c) a metallic Fermionic state coexisting with superfluid phase of Bosons, and d) Mott insulating 
state of Fermions and Bosons. We chart out the phase diagram of the system and provide analytical 
expressions for the phase boundaries within mean-field theory. We demonstrate that the transition 
between these phases are generically first order with the exception of that between the supersolid 
and the Mott states which is a continuous quantum phase transition. We also obtain the low-energy 
collective excitations of the system in these phases. Finally, we study the particle-hole excitations 
in the Mott insulating phase and use it to determine the dynamical critical exponent z for the 
supersolid-Mott insulator transition. We discuss experiments which can test our theory. 

PACS numbers: 67.60.Fp, 64.70.Tg, 73.22.Gk, 73.22.Lp 



I. INTRODUCTION 

Recent experiments on ultracold trapped atomic gases 
have opened a new window onto the phases of quan- 
tum matter—. A gas of Bosonic atoms in an optical or 
magnetic trap has been reversibly tuned between super- 
fluid (SF) and insulating ground states by varying the 
strength of a periodic potential produced by standing 
optical waves 1 . This transition has been explained on 
the basis of the Bose-Hubbard model with on-site re- 
pulsive interactions and hopping between nearest neigh- 
boring sites of the lattice^^. Further, theoretical stud- 
ies of Bosonic atoms with spin and/or pseudospin have 
also been undertaken^. These studies have revealed 
a variety of interesting Motlii and supersolid^ phases 
and superfluid-insulator transitions 3 in these systems. 
On the Fermionic side, the experimental studies have 
mainly concentrated on the observation of paired super- 
fluid states^ and the BCS-BEC crossover in such systems 
near a Feshbach resonance^. More recently, it has been 
possible to generate mixtures of Fermionic and Bosonic 
atoms in a trap^. Several theoretical studies followed 
soon, which established such Bose-Fermi mixtures to be 
interesting physical systems in their own righ t 10 ' 11 , ex- 
hibiting exciting Mott phases in the presence of an opti- 
cal lattice. 

Many of the earlier studies of Bose-Fermi mixtures 
has been restricted to one-dimensional (ID) systems^ or 
have concentrated on regimes where the coupling between 
Bosons and Fermions are weak 13 . The existence of a su- 
persolid (SS) phase in these system in such a weak cou- 
pling regime has also been predicted^. Other work s 10 ' 11 
which have looked at the strong coupling regime have 
restricted themselves to integer filling factors of Bosons 
and Fermions (per spin) and have therefore not addressed 



the phenomenon of translational symmetry breaking and 
possible associated SS phases in the strongly interacting 
regimes of these systems. More recently, however, the au- 
thors of Ref . lla have studied a mixture of spinless softcore 
Bosons with an on-site interaction U and spin-polarized 
non-interacting Fermions at half-filling in a 3D optical 
lattice using dynamical mean- field theory (DMFT) . Sev- 
eral interesting phases, including a SS phase of Bosons 
and charge-density wave(CDW) states of Fermions, have 
been found in Ref. flEI . 

In this work, we study a mixture of hardcore spin- 
less Bosons and spin-polarized interacting Fermions in 
an optical lattice at half-filling using a slave-boson mean- 
field technique. We concentrate on the case where both 
the Bosons and the Fermions have a nearest-neighbor 
density-density repulsive interaction in addition to the 
usual on-site interaction term between them. We pro- 
vide an analytical, albeit mean-field, phase diagram for 
the system and demonstrate that the ground state of such 
a system consists of four distinct phases, namely, a) a 
Mott insulating (MI) phase where both Fermions and 
Bosons are localized at the lattice sites, b) a metal+SF 
phase where the Fermions are in a metallic phase with 
a gapless Fermi surface and the Bosons are in a super- 
fluid state, c) a SS phase where the Bosons are in the 
SS state while the Fermions are localized at the lattice 
site, and d) a CDW+MI phase of coexisting Fermions 
with weak density wave order along with Mott insulat- 
ing Bosons. We show, within mean-field theory, that the 
transition between these phases are generically first or- 
der with the exception of that between the SS and the 
MI phases which is a continuous quantum phase tran- 
sition. We also obtain the low-energy collective modes 
in the metal+SF, CDW+MI and the SS phases and 
demonstrate that they have linear dispersions with defi- 
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nite group velocities. Further, in the MI phase which has 
no gapless modes, we find the dispersion of the gapped 
particle-hole excitations and use it to determine the dy- 
namical critical exponent z for the continuous MI-SS 
transition. We also discuss realistic experiments which 
can test our theory. 

The plan of rest of the paper is as follows. In Sec. 
Ill A[ we develop a slave-Boson mean-field theory for the 
system. This is followed by Sec. Ill B[ where the mean- 
field phase diagram is charted out. In Sec. IIII( we ob- 
tain the low-energy collective modes of the system in the 
metal+SF, CDW+MI and the SS phases. This is followed 
by Sec. lIVI where we discuss the gapped particle-hole ex- 
citations of the MI phase and use it to determine z for 
the SS-MI transition. We discuss relevant experiments 
which can test our theory and conclude in Sec. fVl 



II. SLAVE-BOSON MEAN-FIELD THEORY 
A. Mean-field equations 



[4] can be written as 

H' F = -tpY^dHa^ +V F J2 n >" 

{ij) (y) 

+ (5) 

i 

H' FB = U^nfnf (6) 

i 

where we have implemented the constraint nf = nf using 
Lagrange multipliers Xi at each site, and have used the 
fact that nf — nf — nf < 1 at each site i. We note that 
the Hamiltonian H' = Hb + H' F + H' F is exact and is 
completely equivalent to H (Eq. [1} . 

To make further progress, we proceed with mean-field 
approximation of H' . To this end, we first decompose 
the quartic hopping term in Eq. [5] using 

E + t \ — ^ t \ — A -^^2^1/ \ 

<y> y (ij) 



The Hamiltonian of the Bose-Fermi mixture in a d- 
dimensional hypercubic lattice is described by the Hamil- 
tonian 

H = H F + H B +H FB (1) 

H F = -t F J^cl Cj + V F J^n(nf (2) 

{ij) {ij) 
H B = -tg^blbj + VB^nfnf (3) 

{ij) {ij) 

H FB = Uj^nfnf (4) 

i 

where Cj (bi) and nf = qcj (nf = b\bi) denote the anni- 
hilation and number operators for Fermions (Bosons) at 
site i, V f (Vb) and t F (tB) denote the nearest- neighbor 
interaction strengths and hopping amplitudes for the 
Fermions(Bosons) respectively, U represents the ampli- 
tude for on-site interaction between the Fermions and the 
Bosons, and represents sum over nearest neighbor 

ij pairs on the lattice. In what follows, we shall study 
the Hamiltonian at half-filling. We note at the outset 
that this constraint of half-filling implies nf,nf < 1 at 
each site. 

To obtain an analytical understanding of the phases of 
these model, we first introduce a slave-boson representa- 
tion for the Fermions: Cj = ajdj, where a.; denotes annihi- 
lation operator for the slave-bosons and di represents the 
annihilation operator for pseudo- fermions. We note that 
in this representation, the anticommutation relation for 
the Fermionic operator [cjcj]+ = <5y , where Sij denotes 
the Kronecker delta, enforces the constraint nf = nf on 
each site. In terms of these slave-bosons and pseudo- 
fermions, the Hamiltonians H F and H F b in Eqs. [5] and 



where N is the number of sites in the lattice and the 
hopping amplitudes t\ and ti are given by 

= |E<^i>o- * a = ^£<°I°i>°- (8) 

{ij) ^3) 

Here (..)□ denotes average with respect to the ground 
state of the system. Next, we use a mean-field approxi- 
mation for the constraint term and approximate the La- 
grange multiplier field = Ao + A(— 1)\ where for sake 
of definiteness, we take i = i\ + %2 + ... + id to be even 
for A sublattice sites. Such an ansatz for Ai is motivated 
by the fact that it is the simplest mean-field ansatz that 
preserves the basic symmetries of the problem and, at the 
same time, allows for translational symmetry breaking in 
the pseudo-fermion sector. With these approximations, 
the mean-field Hamiltonian for the system can be written 
as 

Hmf = -i 2 ^4d,+A^(-ir«-nf) 

{ij) 1 

+ E (v F n?n? - tiataj) + J2 Un >i 

{ij) » 

+ Y J {v B nfnf-t B b\b J ) + ^ (9) 

{ij) F 

In writing Eq.[9l we have ignored the term Ao ^2i(nf~nf) 
since it merely renormalizes the chemical potential of the 
fermions and thus behave like a constant as long as we 
restrict ourselves to half-filling. 

To obtain the ground state energy corresponding to 
this mean-field Hamiltonian, we now use a variational 



3 



ansatz for the ground-state wavefunction 



Ha) 
Hb) 
|FS) 



n iv»a> ® n i^) ® i fs ) 



(cos(6»)|n s = 0) + sin(6»)|n s = 1)) 
<g> (cos( 7 )K = 0) + sin( 7 )|n a = 1)) 
(cos(0)\n B = 1) +sin(f9)|n s = 0)) 
® (cos(7)|n a = 1) + sin( 7 )|n a = 0)) 

n^(|k|-M4|0) 



(10) 



where kp denotes the Fermi wavevector for the pseudo- 
fermions and 6 and 7 are variational parameters which 
has to be determined by minimizing the ground state 
energy. Note that the variational wavefunction given by 
Eq. [Til] has a two sublattice structure which allows for the 
possibility of translational symmetry broken phases. For 
the current system, where the Fermions and the Bosons 
both interact via nearest-neighbor density-density inter- 
action terms, Eq. [TU] is the simplest possible mean-field 
variational wavefunction which respects all the symme- 
tries of the Hamiltonian. We would like to point out 
here that capturing the phases of a Bose-Fermi mixture 
with such interaction terms is beyond the scope of any 
single-site mean-field theory including single site DMFT. 

The variational mean-field energy E v — (^\Hmf\^) 
of the system can now be easily obtained and is given by 



Ntp 



-- [Z'sin 2 (27) + Zsin 2 (2(9)] 




:os 2 ( 7 ) cos 2 (A) + sin 2 (7) sin 2 (0)] 
cos(2 7 ) + A(FS|^(-l) i 4d i |FS) 



Z' = {t 1 -V F )/tp, Z={t B -V B )/t F (11) 

where we have used t%/tp = sin 2 (27), A' = A/tF, 
U' = U/tp, and ti/tp has to be determined from Eq. [8] 
The corresponding mean-field equations which determine 
the ground-state values of the variational parameters are 
given by dE v /dd = dE v /dj = dE v /dA' = and yields 



sin(20) 



U' 

— cos(27) + 2Zdcos(26) 



sin(27) 



U' 

— cos(2(9) + 2Z'dcos(2 7 ) - A' 



cos(2 7 ) + -(FS|^(-l) i 4d i |FS) 






(12) 



Next, we evaluate the effective hopping amplitude 

ti and 2/N(FS\Y /i (-iy4^\ FS )- To this end > first > 
let us consider the pseudo-fermion Hamiltonian H = 



AEi(-i) 



The energy spectrum 



of H is given by ±E{k), where E(k) = ^(k) 2 + A 2 , 
e(k) = — 2t2 J2i=i d cos (kia) is the energy dispersion of 



free fermions in a hypercubic lattice in d dimensions, and 
a is the lattice spacing which we shall, from now on, set 
to unity. The density of states (DOS) corresponding to 
these fermions are therefore given by 



P(E) 



Me J 



Pa 



E 



i/||<*>-X><*<>] 



(13) 



where po denotes the DOS of free fermions with tight- 
binding dispersion on a hypercubic lattice, e = e/2^2, 
and we have used the relation p(E)dE = po(e)de. It is 
convenient to use Eq.[lj5]to express the expectation values 
over pseudo-fermion ground states in Eqs. [51 and [T2l and 
one obtains 



cos(27) 



-All (A) 
L°r 1 



h = t 2 h 
TT^A)(S)de 



(14) 



J_°iPo(e)de 



J_iepo(g)rfe 
/"iPo(e)rfe 



h = 



(15) 



Eqs. [12l [14] and [15] denotes the complete set of mean- 
field equations which can be now solved to determine the 
mean-field phase diagram. 



B. Phase diagram 

Eqs. [T2l and [T4l can be easily solved numerically to ob- 
tain the mean-field phase diagram for the system. How- 
ever, before resorting to numerics, we provide a qualita- 
tive discussion of the nature of the phases. 

We find that Eqs. [T^] and [JJ] yields four distinct so- 
lutions which corresponds to four possible phases of the 
system. First, we find a MI phase with broken transla- 
tional symmetry where both the fermions and the Bosons 
are localized. Such a phase corresponds to the solution* 6 



= 0, 7 = tt/2, A 



(16) 



Note that the divergence of A corresponds to ti — > 
which in turn ensures that cos(27) = — 1. Such a MI state 
corresponds to a intertwined checkerboard density-wave 
pattern where the Fermions are localized in sublattice A 



(n? = sin 2 (7) 



= 1 for i € A) and the Bosons are 



localized in sublattice B (nf = cos 2 {6) = 1 for i € B) 
The mean- field energy of this state is E\ =0. 

Second, we find a SS phase, where the Bosons are in 
a supersolid phase with coexisting density-wave and su- 
perfluid order and the Fermions are localized in a Mott 
phase. Such a state corresponds to the solution 



cos(2(9) 



U' 

AZd' 



7 = tt/2, 



00 



(17) 
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Such a state has (6) = sin(20)/2 + 1 and ((-l^&tft,) = 
— cos(20) ^ and thus corresponds to a SS phase for the 
Bosons. Note that the realization of this state necessarily 
requires U'/AZd < 1. For U'/AZd = 1, = and we 
recover the MI state where (b) =0. The energy of the 
SS state is per site given by 



Zdt F ( U' 



\AZd 



- 1 



(18) 



Third, we find the MI+CDW state where the fermions 
show weak density-wave oscillations whereas the Bosons 
are localized in the MI state. This corresponds to the 
solution 



0, 7-70^0,^/2 



(19) 



where 70 and A are to be determined from a numerical 
solution of the mean-field equations 



cos(27o) 



U' 
AZ'd 



(1-2A/J7) = -A7i(A). (20) 



The energy of this state per site is given by 
U cos 2 (70) 



£?a — 



/ AZ'd . A 

^l-_sin 2 ( 7 o)J (21) 



Finally, we find the state in which the superfluid 
Bosons coexist with metallic Fermions. This corresponds 
to the solution 



7 = 7r/4 A = 



(22) 



Note that such a phase has (b) — sin(2#)/2 = 1 and 
((— iyb\bi) = — 008(26*) = so that the Bosons are in an 
uniform superfluid state. Also, A = and ((— l) l a\ai) — 
— cos(27) = in this state indicating that the Fermions 
are in a gapless uniform metallic state. The energy of 
this state per site is given by 



(23) 



The phase boundaries corresponding to these phases 
can be analytically computed using Eqs. [TH] [21] and [23l 
provided 70 and t\ (which determines Z') are obtained 
from numerical solutions of Eqs. [501 and Eq. [5] For the 
MI phase to occur, we must have E2, E3, £4 > E\ = 
which yields the conditions 



2d 



U' 



1--(Z + Z'))>0, Z<0,nd—^<1 



A\Z\d 



AZ'd 



U' 



-sin 2 ( 7o )<l (24) 



Note that the condition U' / (A\Z\d) < 1 which is nec- 
essary for the realization of the SS phase has to be si- 
multaneously satisfied with the condition Z < to make 



sure that the SS phase is actually a competing candidate 
to the MI state. The MI phase can indeed be realized 
in the parameter regime Z > provided U' > A\Z\d. 
The first condition (Z + Z') < U'/2d shows that the MI 
phase is favored over the metal+SF phase for large U/d 
and predicts a linear phase boundary in the U' — Z plane 
U' = 2d(Z + Z') with a slope of 2d and intercept of 2dZ' 
between these two phases. Note that the MI phase always 
wins over the metal+SF phase if the nearest-neighbor 
interactions between the Bosons and Fermions are large 
compared to their hopping amplitudes making Z + Z' 
negative. The final condition sin 2 (70) < 1 indicates 
that the phase boundary between the MI and MI+CDW 
phases is independent of Z. The former phase is favored 
over the latter for larger U and smaller Z' . 

Similarly for the SS phase to occur one needs 
U'/{A\Z\d) < 1 and E 2 < £4, £3, £4, which yields 



Zd 



« (AZ'd n . . \ / U' 

cos 2 ^— sm 2 (7o)-l) <(^-l 



Z > 0, 



AZ'd 



U 



— sin (70) < 1, 



U' 



A^fZZ'd 



> 1 



(25) 



We note that the SS phase is favored when the nearest- 
neighbor interaction between the Bosons are weak com- 
pared to their hopping amplitudes making Z positive 
and when U' is small enough so that U' / (A\Z\d) < 1. 
Also, from the conditions in Eq. [55] (obtained using 
E 2 < £4, Ei), we note that for a given Z' and d, the 
boundary between the metal+SF and the SS phases is a 
parabola in the U' — Z plane given by U 2 = 16ZZ'd 2 
while that between the SS and the MI state is a line given 
by U' = AZd. 

Finally, the condition for occurrence of the metal+SF 
phase is given by £4 < Ei,E%, E3 and is given by 



2d 



l-—(Z + Z')) <0, 



U' 



A^JZZ'd 



( 2d, .A n, J AZ'd 2 , . 

fl-— (Z + Z')j <2cos 2 ( 7 o) (^l--p--sin 2 ( 70 ) 



(26) 



The last condition in Eq. l26l determines the phase bound- 
ary between the metal+SF and the MI+CDW phases 
which depends on value of 70. However, numerically, we 
find that for U ~ 70 — 7r/4, and in this regime, the 
phase boundary between these phases occurs at Z ~ for 
all Z' and d. Note that strictly at U = 0, the Fermionic 
state is metallic; however a CDW gap opens up in the 
Fermionic spectrum for an infinitesimal finite U' . 

To verify the above-mentioned qualitative arguments 
and to find a precise phase diagram for the system, we 
numerically solve Eqs. [12] fT4l and [T51 for d = 2 and for 
representative values Vp/t F = 0, 0.5. We plot the ground 
state phase diagram as a function of Z and U' in Figs. 
Q] and [U We find that the numerical results agree well 
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FIG. 1: Ground state phase diagram as a function of Z and U' 
for noninteracting Fermions (Vf = 0). The phase boundaries 
coincides with the analytical mean-field phase boundaries (see 
text for details). 
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FIG. 2: Same as in Fig. Q]but for V F /t F = 0.5. The interac- 
tion between the Fermions favors the SS phase as can be seen 
by comparing Fig. [T] and [2] 



with the qualitative arguments. Figs. [T] and indicate 
that the phase boundary between the CDW+MI and MI 
phases is independent of Z as noted earlier. The linear 
and the parabolic nature of the phase boundaries between 
the MI and SS phases and the SS and metal+SF phases 
respectively can also be easily verified from the Figs, and 
are in accordance with the qualitative discussion. We 
note that one of the effects of nearest-neighbor repulsion 
between the Fermions is to enhance the SS phase which 
occupies a larger region of phase space in Fig. [2] (Vf ftp = 
0.5) than in Fig. Q] (Vf = 0). Such an interaction, for 
Z < 0, also favors the Mott phase over the CDW+MI 
phase as can also be seen from Figs. [I] and [2] 

To determine the nature of transition between the dif- 
ferent phases, we plot the ground state values 8 as a 
function of Z for U' = 10 and V F /t F = 0.5 in Fig. H 
Such a plot clearly shows that the transition between the 
metal+SF and the SS phases is, within the mean-field 
theory considered here, first order and is accompanied by 
a jump in the value of 9. In contrast, the SS-MI transition 
turns out to be continuous. A similar plot of ground state 
values 7 as a function of U' for Z/tp = —2 and Vf = , 
shown in Fig.[H indicates that the transition between the 
CDW+MI and the MI phase is also discontinuous and is 
accompanied by a jump in the ground state value of 7. 

Next, we compare our phase diagram with that ob- 
tained from DMFT in Ref. 15. This can be done in 
the regime of large positive Z (which correspond to 
VbI^b — > 0) which was the case treated in Ref. [H. We 
find that the two phase diagrams qualitatively agree in 
the sense that both yield SS and MI phases in these limit. 
The difference lies in the fact that our mean-field predicts 
a second-order transition between the two phases whereas 



DMFT yields a narrow region of coexistence. This is pre- 
sumably an effect of quantum fluctuation which is not 
captured within the mean-field theory. In addition, we 
also find a region of metal+SF phase at low U which was 
not seen in Ref. Il5l . 

Finally, we would like to point out that the slave-boson 
mean-field phase diagram obtained above yields qualita- 
tively correct phase diagram, but not a quantitatively 
correct one. This can be most clearly seen by noting that 
our Hamiltonian reduces to an effective Falicov-Kimball 
(FK) modelil in the limit Z = t B = V B = Vf = 0. This 
is most easily seen by writing our starting Hamiltonian 
H (Eq. H} for t B = V B == V F = 

H FK = -t F J2 c lcj + Uj2chnf (27) 

(ij) 1 

At half-filling the Bosons are localized in the B sublattice 
so that nf = (1 — (— l) l )/2. Thus the Fermions have a 
CDW instability even for an infinitesimal U due to nest- 
ing for half-filling on a square lattice. Consequently, the 
FK model at half-filling is insulating for any infinitesi- 
mal U, as known from several earlier studies^. Such a 
CDW instability, which can be easily captured by weak- 
coupling mean-field theory, is not straightforward to ob- 
tain in our strong coupling slave-boson mean-field ap- 
proach which predicts a finite critical U for the transition 
from metal+SF to the MI phase. Note however that the 
slave-boson mean-field theory does predict a CDW+MI 
state for weak U, but for small negative Z, as can be seen 
from Fig. [TJ This indicates that the phase diagram ob- 
tained has qualitatively, but not quantitatively, correct 
features. 
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FIG. 3: (Color online) Variation of ground state value of 6 
with Z for U' — 10 and Vf jtv = 0.5. The discontinuity in 
9 at Z = 3.38 indicates a first order transition between the 
metal+SF and the SS phases. The transition between the SS 
and the MI phases is continuous. 




FIG. 4: (Color online) Variation of the ground state value 
of 7 with U' for Z — —2. and Vf — 0. The discontinuity 
in 7 at U' — 2 indicates a first order transition between the 
CDW+MI and the MI phases. 



III. COLLECTIVE MODES 



The gapped modes corresponds to particle-hole excita- 
tions in the CDW+MI, SS and the Mott states. In this 
section, we concentrate on the gapless collective modes 
in the CDW+MI, SS, and the metal+SF phases. 

To obtain the dispersion, we first consider a time- 
dependent variational wave-function 

\9 d (t)) = (cos(0i)|nf =O)+sin(0 i )e- i * i |nf = 1)) 

® (cos(7i)|nf = 0) +sin( 7i )e-*|nf = 1)) 

(28) 

where 0j, xu Ti, and are space-time dependent fields. 
Note that in the static limit, the ground-state of the sys- 
tem corresponds to 9t = 9io = 9q(it/2 — 9q), 7, = 7,0 = 
7o(7r/2 — 70) for i G A(B) sites, and Xi — ^i — 0, so that 
|*d(*)) reduces to |*). 

The Lagrangian C = *£i(*d(t)\id t -H' + » B £. nf + 
Mf J2i n f \^d) can now be computed using the variational 
wave-function and one obtains 



C = 53[^Xism 2 (0i)+$&sm 2 (7i) 

i 

-C/sin 2 (^)sin 2 ( 7i ) 

+ S [4 sin ( 2 ^) sin(2^) cos( Xl - Xj) 

(ij) 

+ii sin(27i) sin(27_,-) cos(0 4 - 4> 3 ) 

-{V F sin 2 ( 74 ) sin 2 (7,) + V B sin 2 ^) sin 2 (0,-)}) 

+ ^[A{-sin 2 ( 7l )(-l) l + «)} 



sin 2 (6'i) + hf sin 2 (7i) 



(29) 



where 



) = (FS\(-l) i 4d i \FS)/N is the Fermion 
number density different on A and B sublattices and all 
time dependence of the fields are kept implicit for the 
sake of clarity. To obtain the collective modes, we now 
write 9i{t) = 9 l0 + £0,(t), 7,(t) = % + #7i(*)i and ex- 
pand the Lagrangian to quadratic order in <J0j(i), (i), 
</>i(i) and Xi(t)- Then a variation of this Lagrangian with 
respect to S9i(t), #7i(i), </>i(i) and an d consequent 

adjustment of values of the parameters /is and fj,p fol- 
lowing Ref. El yields the equations for the low-energy 
collective modes (we set h = 1 from now on) 



The phases of the Bose-Fermi mixture discussed in ^7k + tidsin(2j ) (1 - c(k)) <^k = (30) 

the previous section allows for two types of excitations. <9t</>k — ai(k)^7k — C/sin(20 o )<50k = (31) 

The first type is the low-energy gapless collective modes d t S9 k + t B dam(29 ) (1 - c(k)) X k = (32) 

that are present in the metal+SF, CDW+MI, and the SS * J ' ) , ' 

phases of the system. These are the gapless Goldstone dtXk " a 2(k)^ k - t/sm(2 7o )<h k = (33) 
modes corresponding to the Boson and the slave-Boson 

fields. For all of these modes, the pseudo-Fermion sector where we have taken the Fourier transform of all the 

remain gapped and do not contribute to their dispersion. fields, c(k) = Ylj=i d cos(kj)/d, and ct\ and a-i are given 
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by 



«i(k) 



a 2 (k) 



+ C (k)(^sin 2 (2 7o ) 



sin(27o) 
+ cos 2 (2 7o ) 



4 ' fld r i + c(k)(^sin 2 (20 o ) 



sin(2fl 
+ cos 2 (26» ) 



(34) 



It is important to note that Eqs. 13171 and I3T1 holds when 
7o ^ 0,tt/2 while Eqs. EH and [33] holds when 9 ^ 0,tt/2. 
Thus none of these equations are valid in the MI phase 
which do not support any low-energy collective modes. 
The gapped modes of the MI phases will be obtained in 
the next section. 

In the SS phase where cos(26» ) = U'/(AZd) £ 0, 1 
and 7o = tt/2, the collective mode corresponds to the 
low-energy excitations of the Bosons and are given by 
Eqs. 1321 and [331 A simple set of standard manipulations 
of these equations yield the dispersion of the collective 
modes lo 2 = 2u 2 (k)(l — c(k)), where 



vi(k) = t B dsin{29 )a 2 (k)/2 



(35) 



Note that for low momentum, we get a gapless linearly 
dispersing collective mode with velocity v ss = Vi(k = 0). 
Similarly for the MI+CDW phase, where 9 = and 
7 = 70, the collective mode corresponds to the low energy 
excitations of the pseudo-bosons and can be obtained by 
solving Eqs. 1301 and 1311 Since the pseudo-Fermion sector 
is always gapped in this phase (A ^ 0), the collective 
mode here corresponds to the density-wave mode of the 
real Fermions. These modes have linear dispersion uj 2 — 
2uf(k)(l - c(k)) where 



vl(k) = iidsin(27o)ai(k)/2 



(36) 



Thus for low momenta, we again get a gapless lin- 
early dispersing collective mode with velocity vcdw = 
«a(k = 0). 

Finally for the metal+SS phase, all the Eqs. [301 [331 
hold and they need to be solved simultaneously. In 
this phase since 70 = 9q = tt/4, we find that ai(k) = 
4M[1 + c(k)V F /ti] and a 2 (k) = At B d[l + c{k)V B /t B ]. 
Solving these equations, one finds two collective modes 



with linear dispersions lo\ 
are given by 



2w|(k(l-c(k)) where u±(k) 



w ±(k) = - [(a 1 (k)t 1 d + a 2 (k)t B d) 



±y/(ai(k)hd - a 2 {k)t B d) 2 + 16(?7i s tid) s 



(37) 



These collective modes result from the hybridization of 
the Bogoliubov modes of the Bosons and the density- 
wave modes for the metallic Fermions. This fact can be 
easily checked by putting U = in Eq. [37] by which 
one can retrieve these modes with velocities v B (k) = 



a 2 (k)t B d/2 for the Bosons and v F (k) = ax(k)txd/2 for 
the Fermions. As U' increases, the hybridization be- 
tween these modes become stronger until the velocity 
V-(k = 7r) touches zero at U' = Ay/\Z'Z\d which is pre- 
cisely the condition for the metal+SF phase to become 
unstable to the SS phase. 



IV. GAPPED MODES IN THE MI PHASE 

The MI phase, in contrast to the other three phases of 
the system, do not support a gapless mode. The lowest- 
lying excitations of such a state with conserved number 
density arc particle-hole excitations. Such excitations, 
are of two types. The first type, shown in second panel 
of Fig. [5] involves particle and hole excitations that spa- 
tially well-separated while the second type, shown in sec- 
ond panel of Fig. [6] involves particle and hole excitations 
in nearest-neighbor sites which forms a dipole. In what 
follows, we first compute the energies of both these exci- 
tations using perturbation theory up to second order in 
ifl if /Vb if which are supposed to small in the MI phase. 

Such an energy estimate can be easily carried out by 
strong-coupling perturbation theory developed in Ref. [2(J 
in context single species Bose-Hubbard model. The gen- 
eralization is largely trivial, except for one important de- 
tail. In the standard Bose-Hubbard model studied in 
Ref. any particle/hole excitation could have lowering 
of energy via nearest-neighbor hopping which is 0(t/U) 
process. In contrast, as can be seen from Figs. [S] and 
[6] it is not possible for the particle-hole or dipole ex- 
citations to directly hop to the next site since such a 
direct hop always take us out the low energy manifold 
of states in the Mott limit. In particular, we note that 
any kinetic energy gain of the particle-hole or dipole ex- 
citation must occur via hopping of the partice/hole to 
the second-neighbor sites and hence necessarily leads to 
0(t 2 /V 2 ) energy gain. 

We first compute the excitation energy of the 
Bosonic(Fermionic) particle-hole pair when they are far 

apart. The on-site energy of creating such a pair is 

B ip 

^on-site = ^dV B / F + U while the energy-lowering due 
to hopping of each of the particle and the hole is given 
^ ^pping = -2d(2d-l)t 2 B/F /{2(2d~2)V B/F + U}. The 
energy of the Mott state to second order in perturbation 
theory is E M i = -d{t 2 B /[2{2d - 1)V B + U] + t 2 F /[2(2d - 
1)Vf + U]) so that the excitation energy of the particle- 
hole pair is 



E 



5/F 
p—h 



AdV, 



B/F 



U 



dt% 



M(2d-l)t 2 B/F 
2{2d-2)V B/F + U 
dt 2 F 



7 (38) 



2{2d-l)V B + U 2(2d-l)V F + U y 
We note that in the limit of large d, where the mean-field 
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FIG. 5: (Color online) Cartoon representation of the MI state. 
A) The Mott state at half-filling. The red filled circles repre- 
sent Bosons and the empty blue circles indicate Fermions. B) 
A particle and a hole excitation which are not nearest neigh- 
bors. C) An intermediate virtual high energy state which as- 
sists hopping of holes. D) A state where the hole has hopped 
to the next-neighbor site. This state is identical in energy to 
the state B. 



tn 



o- 



-e- 



(A) 



(C) 



t 



-Q 



B53 

(B) 

tn 



o 



e-m- 

(D) 



FIG. 6: (Color online) Cartoon representation of the MI state 
and the associated dipole excitations. All symbols are same 
as in Fig. [5] A) The Mott state at half-filling. B) A dipole 
excitation over the Mott state C) An intermediate virtual 
high energy state which assists hopping of dipoles. D) A state 
where the dipole has hopped to an adjacent link. This state 
is identical in energy to the state B when Vf =Vb- 



results are expected to be accurate, we have 



E 



B/F 



AdV } 



B/F 



u ■ 



8d 2 t 2 B/F 



B/F 



u 



(39) 



The Mott state is destabilized in favor of the SS phase 

when E B p Ll = 0. 

Next, we compute the excitation energy of the 
Bosonic/Fermionic dipole state. We are going to do this 
in the limit of Vp = Vb = V. We note at the outset that 
once such a dipolar excitation is created, it remains sta- 
ble, i.e., the hole can not hop away arbitrarily far away 
from the particle. It can be easily verified from Fig. [6] 
that such hoppings generate higher-energy end states and 
takes one out of the low energy manifold of states. 

The on-site energy cost for creating such an excitation 
is E on - s ite — 2(2d — l)V + U while the hopping process, 
shown schematically in Fig. [SJ necessarily involves hop- 
ping of both Fermions and Bosons and leads to an energy 
gain of E£ a?piag = t B t F /2(2d- l)V . Thus the net energy 
of such a dipole excitation is given by 



E, 



dipole 



2 (2d - 1)V + U 



t B tp 



dt% 



2{2d- 1)V 
dtp 



2(2d-l)V + U 2(2d-l)V + U 



(40) 



We note that for d ^> 1, where our mean- field theory 
holds, it is always energetically favorable to create parti- 
cles and holes well-separated since E^ p F ping < E$ ing . 
However, the dipolar excitations may becomes favor- 
able in low dimension and for large U/V. In this case, 



^SJing » popping for U » V and in this limit, the 
dipolar excitations would be preferred in destabilizing the 
MI phase. We shall not discuss this issue here any further 
since this is clearly beyond the scope of our mean-field 
theory. 

Finally, we compute the dispersion of the gapped 
particle-hole excitations within mean-field theory where 
the particle and the hole are spatially well-separated and 
do not interact. To this end, we temporarily relax the 
constraint of conservation of particle-number and con- 
sider the energy of excitations of adding a particle E p 
and a hole Eh to the Mott state. The physical particle- 
hole excitation energy can then be computed from E p h = 
E p + Eh- To compute the energy of these particel/hole 
excitations, we adapt a time-dependent variational ap- 
proach as done in Ref. [HI for single species Bosons in 
an optical lattice. We begin with the variational wave- 
function of Bosons and slave Bosons(Fermions) given by 



|V(t)> = x I'M*)) 

\Mt)) - / o Q (i)|n s = 0) + /f(t)|n s = l) 

= 9S(t)\n a - 0) + g?(t)\n a = 1) (41) 

where a = A, B denotes sublattice indices. The co- 
efficients / and g satisfy the normalization condition 
l/o| 2 + |/i| 2 = 1 and M 2 + \ 9l \ 2 - 1. We note that 



at equilibrium /q = sin 9 = 0, ff = 



cos 7 



0, g? 



sni7 



1. 



fB 
J 0(1) 



COS ( 

= ft 



(0) 



= 1, 

, and 



00(1) = #i(o) for the MI Phase. 

The Lagrangian of the Bose-Fermi mixture can then 
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be written as 

L ' = E i [fojfoj + fij.fi] + 9oj90j + g*j9ij 

3 

- E [ - ^fufoifojhj ~ h9ti9oi9tj9ij 

(ij) 

+V B \fu\ 2 \hj\ 2 + V F \gu\ 2 \gi 1 \ 2 

-tf£i/iii 2 M 2 +w£i/«i 

i i 
i i 

-£^(l<7io| 2 + |sii| 2 -l 



/Oil 



I/hI 2 



1) 



(42) 



where Xi and i/j are variational parameters used for im- 
plementing the constraint whose values are to be deter- 
mined from proper choice of the saddle point which in 
the MI phase yields A^ = fib and Xb = for the Bosons 
and va = and vb = for the slave bosons. 

The saddle-point equations for the variational coeffi- 
cients fi(t) and gi(t) then reads 







ifoi + tsfu £ fijfoj 

(i> 

igoi + hgu E 9*j90j - "igoi = 



(?) 



<B/oiE/oi/ii+2Vfl/nEl/bl 



Xifii 



U) 



W F gii E 



+^9ii|/ii| 2 - M/Sli] - ftfu = (43) 



Next we implement the two sublattice structure, shift 



to momentum and frequency space, and expand / ak 



A/B 



j-A/B , A/B r A/B 

/ a and 5ak = <S 3ak 



-ga^ B where a = 0, 1. 



Note that since 5f A / B and 5g A / B corresponds to devi- 
ation of particle number of the MI state, these disper- 
sion corresponding to their eigenmodes must represent 
the particle and hole excitations over the MI phase. In 
the MI phase, we find that the equation of motions for 
the Bosons and the pseudobosons decouple at linear or- 
der 5f^ B and 5g~^ B . For the Bosons, we obtain, to 
linear order in Sf^ B 



-2dt B c(k)Sf* B + X A Sf A 
-2t B dc(k)5f& + 2V B zSfg 
+U6f* k - ( M6 - X B )6f B (44) 



sponding to particle and hole excitations 

U 



E 



■ P (h) 



-(-)(2dV B + -- M 6) 



(2dV B + j) 2 - (2t B dc(k)) 2 (45) 



The energy of a particle-hole excitation which conserves 
particle number is therefore obtained by adding E p and 
Eh and is given by 



E v _ h = 2\[2dV 1 



|)2 - (2t B dc(k)Y 



(46) 



which yields two physical excitation dispersion corre- 



Note that E p _h vanishes along the line Z = U' /Ad which 
agrees with the mean-field result for the SS-MI phase 
boundary. Also, expanding Eq. US] to 0(t B ) for k = 
leads exactly to Eq. [39] which shows that the second- 
order perturbation theory discussed earlier agrees to the 
present calculation in the high d limit. Further, at small 
wave-vector, we find E p _h ~ |k| which shows that the 
SS-MI quantum phase transition has a dynamical critical 
exponent z = 1. Similar dispersion can be obtained for 
the pseudo-bosons by considering collective modes corre- 
sponding to Sg a k. These modes have the same dispersion 
as Eas.HBI and [4*61 with Vb and ts replaced by Vf and t\ 
respectively. 



V. DISCUSSION 

Experimental realization of Bose-Fermi mixtures have 
long been achieved in ultracold atomic systems. These 
mixtures can be easily tuned to a regime where the on- 
site intra-species interaction between both the Fermions 
and the Bosons are large so that they effectively be- 
have as hard-core particles. Thus experimental realiza- 
tion of a Bose-Fermi mixture with Vb = Vf = is rela- 
tively straightforward. However, most such mixtures do 
not have sufficiently large nearest-neighbor repulsion and 
thus it might be difficult to realize mixtures which has 
large Vf or Vb- Some progress in this direction has re- 
cently been made in Ref. |21| . Also, use of spin-polarized 
52 Cr atoms for the Fermionic part of the mixture may 
help since these atoms have significant dipole moment 
which may provide the requisite interaction. 

Once such a Bose-Fermi system is realized, several pre- 
dictions of the present work can be verified by realizable 
experiments that are commonly used for ultracold sys- 
tems. First, we note that since the Bosons are spinless 
and the Fermions are spin-polarized, the Bosonic and the 
Fermionic part of the mixture can be separated by ap- 
plying a standard Stern-Gerlach field during a standard 
time-of-flight experiment as done earlier in Ref. |22] in 
the context of spinor Bosons in optical traps. Such a 
procedure allows us to separately study the momentum 
distribution functions of the Bosons and the Fermions 
using time-of-flight experiments 1 . For the Bosonic cloud, 
the distinction between the SF and the MI phases can be 
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easily done by measuring the presence or absence of co- 
herence peaks in its momentum distribution as measured 
in a standard time-of-flight experiment. The precise na- 
ture of the broken translational symmetry in the MI and 
the SS phases for the Bosons can also be determined by 
studying noise-correlations of the expanding clouds as al- 
ready proposed in Ref.|M Thus, the MI, SS and the SF 
phases for the Bosons can be qualitatively distinguished 
by these experiments. As for the Fermions, the pres- 
ence/absence of a Fermi surface for the Fermions in a 
trap can be easily distinguished in time-of-flight mea- 
surements as performed for ultracold Fermions in Ref. 
l24l . Thus, these experiments should allow one to qualita- 
tively distinguish between all four predicted phases. One 
of the central predictions of our theory is that for any 
finite U, the metallic state of the Fermions shall always 
be accompanied by a SF phase of the Bosons. In terms 
of time-of-flight experiments this means that any mea- 
surement on Fermions which sees a gapless Fermi surface 
shall always be accompanied by corresponding coherence 
peak (and no density wave ordering) for the Bosons. The 
collective modes computed in this work can also be ver- 
ified experimentally using standard inelastic light scat- 
tering experiments^. Such experiments can detect low- 
energy collective modes and should thus detect either two 
( metal+SF phase) or one (SS or MI+CDW phases) lin- 
early dispersing collective mode(s). The MI phase will 
be characterized by absence of any low energy collective 
modes of the system. 

There are several possible extension of our analysis. 



The first and the simplest extension would be to study 
the phases of the Bose-Fermi mixture away from half- 
filling. This would require a more careful handling of 
the chemical potential \ib and [If of the Bosons and the 
Fermions. In particular one would need to determine t\ 
in a self-consistent manner as a function of [ip. Second, it 
would be interesting to look at the phase diagram by re- 
laxing the hard-core constraint on the Bosons by putting 
a finite on-site repulsion between them. Of particular in- 
terest in this respect is to check if the slave-boson mean- 
field can provide any indication of the phase separation 
found in such a system in Refs. [Tlflll . Finally, it would 
be useful to study the phase diagram of mixture of spin- 
polarized Fermions with spin-one and spin-two Bosons 
with nearest-neighbor interactions. Such system are ex- 
pected to have a much richer phase diagram and have 
not been theoretically studied so far. 

To conclude, in this work, we have carried out a slave- 
boson mean-field analysis of a mixture of hardcore spin- 
less Bosons and spin-polarized Fermions in an optical lat- 
tice. Our analysis provides the mean-field phase diagram 
of the system and shows the presence of four distinct 
phases. We have also computed the low-energy collective 
modes of three of these phases (metal+SF, CDW+MI 
and SS) and studied the gapped particle-hole excitation 
of the fourth (MI). We have discussed experiments which 
can be used to test our theory and possible extension of 
our theory to other systems. 

We thank Jim Freericks for drawing our attention to 
the Falicov-Kimball limit of the present model and for 
several useful discussions. 
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